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Abstract.  Self-assembled  3D-islands  formation  during  heteroepitaxial  growth  was  investigated  by 
Monte  Carlo  simulation.  Shwoebel  barriers  for  explanation  of  3D  growth  kinetic  were  suggested. 
Island  size  equalization  during  growth  process  was  observed.  Necessity  of  atomic  flux  from  island 
edge  to  the  upper  layers  for  3D  island  fonnation  was  demonstrated. 


Heteroepitaxial  growth  have  been  utilized  for  the  generation  of  strongly  strained  coherent 
3D  islands,  named  quantum  dots  (QDs)  [1,  2],  The  elastic  energy  of  the  initially  planar 
strained  layer  increases  during  heteroepitaxial  growth.  The  system  can  lower  its  free  energy 
by  atom  transitions  from  the  island  edge  to  the  upper  layer  as  that  leads  to  decrease  of 
contact  area  between  substrate  and  new  layer  [1,  3-6].  Kinetic  of  such  process  could 
be  described  taking  into  consideration  diffusion  peculiarities  when  an  atom  crosses  island 
edges.  Additional  energy  barriers  for  atoms  crossing  monoatomic  step  (Schwoebel  barriers) 
are  suggested  to  be  responsible  for  QD  formation  [7,  8],  Kinetic  problems  of  epitaxial 
growth  with  Shwoebel  barriers  were  considered  theoretically  in  [9],  Increase  of  nucleation 
rate  of  islands  in  the  second  monolayer  was  associated  in  [ 1 0]  with  asymmetry  of  Schwoebel 
barriers  in  SiGe  system. 

Introduction  of  additional  energy  barriers  Es,  for  atom  hops  over  monoatomic  steps  in 
the  model  enable  us  to  investigate  Schwoebel  barriers  influence  on  surface  relief  fonnation 
during  epitaxy  [11].  Diffusion  hop  probability  to  cross  the  step  is  changed  by  the  factor 
P — exp(—Est  / kT)  as  compared  with  diffusion  hops  at  the  same  atomic  level.  Factor  Pup 
changes  probability  of  atom  hops  to  the  upper  layers,  and  Pdown  — to  the  lower  ones.  Phase 
diagram  in  coordinates  Pup  — Pdown  is  presented  in  Fig.  1(a).  Shaded  area  — Stransky- 
Krastanov  growth  mode  - separates  layer-by-layer  regime  (2D)  from  3D-island  growth 


Fig.  1.  (a)  Phase  diagram  indicating  conditions  for  3D  islands  fonnation,  shaded  area  — Stransky- 
Krastanov  growth  mode,  1 — line  Pllp  = Pdown,  2 — X = 40  a.s.,  3 — X = 20  a.s.;  (b)  simulated 
surface  after  1.5  ML  deposition,  X = 20  a.s.,  x = 10. 
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Fig.2.  Number  ofatoms  in  three  Fig.  3.  (a)  CNS  dependence  on  x;  (b)  expectation  time  te 

successive  atomic  layers  of  3D  for  appearance  of  the  nucleus  of  the  second  monolaeyr  on  2D 
island  formed  at  X = 40  a.s.  islands  versus  x ■ 
versus  time  for  x = 2 (fig- 
ure at  the  curve  corresponds  to 
atomic  layer  number),  i * — crit- 
ical (2D)  nucleus  size  (CNS). 


□ 1st  island  ■ 2nd  island 


Fig.  4.  Islands  self-organization  during  growth  process,  X = 40  a.s.:  (a)  histogram  of  size 
redistribution  of  two  islands  with  time,  left  columns  correspond  to  the  moments  of  the  second  layer 
nucleus  creation  on  2D  islands;  top  view  of  model  surface  (200  x 250  a.s.)  after  deposition  0.2  ML 
(b)  and  0.45  ML  (c),  to  the  right  of  each  surface  — its  profile  (b-c). 


mode.  Shaded  area  is  parallel  to  the  line  Fup  = Pdown  which  is  to  say  that  the  change-over 
from  2D  to  3D  mode  (2D  — *■  3D)  is  determined  only  by  the  parameter  x = F’up/F’down  and 
not  by  the  absolute  values  Pup,  Pdown-  Notice  that  for  3D  island  formation  without  wetting 
layer  condition  / > 1 is  necessary.  Increasing  of  migration  adatom  length  X ( D / V),  (I)  — 
surface  diffusion  coefficient,  V — deposition  rate)  leads  to  parallel  shift  of  transition  2D 
— + 3D  to  the  lower  values  Pup,  under  condition  x > 1 • Figure  1(b)  represents  the  model 
surface  with  3D  islands. 

Number  of  atoms  in  three  successive  atomic  layers  of  3D  island  formed  on  ( 1 00)  surface 
at  X = 40  atomic  sites  (a.s.)  versus  time  is  shown  in  Fig.  2.  Fig.  2 demonstrates  stable  20- 
island  formation  till  3D  island  growth  starts.  Size  of  this  2D  island  fluctuates  about  average 
value  i* , that  is  equal  to  critical  nucleus  size  (CNS).  Large  CNS  corresponds  to  active 
exchange  between  islands  and  equilibrium  adatom  gas  due  to  small  binding  energy  between 
neighboring  atoms.  For  equilibrium  conditions  following  relationship  was  obtained  i*  ~ 
( D/V )’/5  that  means  slight  increase  of  i*  with  X increase.  As  one  can  see  in  Fig.  2 CNS 
decreases  in  each  next  atomic  layer. 

The  critical  size  of  island  icr  when  nucleation  of  new  layer  starts  is  of  interest.  Depen- 
dence icr(x)  was  investigated.  Dependencies  i*  and  icr  on  x for  two  X are  presented  in 
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Fig.  5.  Nup/Ndovm  (ratio  of  absolute  number  of  hops  to  the  upper  layers  to  hops  to  the  lower 
layers)  dependence  on  time;  X = 40  a.s.,  1 — x = 2,  2 — x = 10. 

Fig.  3(a).  As  / increases  these  sizes  decrease  approaching  1 for  / — > oc.  From  some  x, 
i*  = icr}  that  means  3D  nucleation  from  the  initial  phase  of  growth.  Expectation  time  te  for 
appearance  of  the  nucleus  of  the  second  monolaeyr  on  2D  islands  versus  x is  represented 
in  Fig.  3(b).  This  dependence  is  exponential  one:  te(x)  = T()  -exp(— or  • / )■  For  A = 40  a.s. 
— to  = 0.017  • to,  a — 0.8,  and  for  X = 20  a.s.  — to  = 0.035  • to,  a = 0.5,  where  to  - 
time  necessary  of  one  monolayer  deposition. 

3D  islands  without  wetting  layer  demonstrate  rather  sharp  size  distribution.  That  is  due 
not  only  to  the  most  synchronous  2D  island  nucleation  in  the  first  atomic  layer.  Lateral 
size  equalization  when  3D  growth  process  starts  (case  of  large  size  fluctuation  of  initial 
2D-islands)  could  be  seen  in  Fig.  4.  Figure  5 demonstrate  increase  with  time  of  hops  up 
to  hops  down  ratio  Nup/Ndov/n.  Such  (Vup/Miown  variation  is  due  to  decrease  of  atoms 
capable  to  attach  into  descending  island  boundary  because  of  their  attachment  to  the  island 
of  upper  layer. 
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